Self-interacting holographic dark energy 
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We investigate a spatially flat Friedmann-Robertson- Walker (FRW) universe where dark matter 
exchanges energy with a self- interacting holographic dark energy (SIHDE). Using the x^-statistical 
method on the Hubble function, we obtain a critical redshift that seems to be consistent with 
both BAO and CMB data. We calculate the theoretical distance modulus for confronting with the 
observational data of SNe la for small redshift z < 0.1 and large redshift 0.1 < z < 1.5. The 
model gets accelerate faster than the ACDM one and it can be a good candidate to alleviate the 
coincidence problem. We also examine the age crisis at high redshift associated with the old quasar 
APM 08279+5255. 
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I. INTRODUCTION 



As it is well known our universe is currently undergo- 
ing an accelerated expansion phase driven by a mysteri- 
ous fuel called dark energy which exerts a negative pres- 
sure tending to drive clusters of galaxies apart. The lat- 
ter fact has been corroborated by many different probes, 
for example the observation of type la supernovae 1- J], 
the data of the large scale structure from SDSS 5143' 
and measurements of the cosmic microwave background 
anisotropy [8l-[l0|. The simplest candidate for the dark 
energy component is a positive cosmological constant A 
jll| - [l3| . Although the prediction of the cold dark mat- 
ter plus cosmological constant (ACDM) model is mostly 
consistent with observational data, the cosmological con- 
stant proposal suffers from at least two puzzles (l3 - l20| . 
The first issue is known as the fine-tuning problem, that 
is, the theoretical prediction of the cosmological constant 
that is obtained as the expectation value of quantum 
fields differs from its cosmic observed value by 120 or- 
ders of magnitude. The measured cosmological constant 
in our universe is tiny but not zero, and if it were much 
larger, galaxies could not have formed [l^l- The second 
point of debate concerns the cosmic coincidence problem: 
why we observe that the fractional densities of dark mat- 
ter and cosmological constant are about the same order 
of magnitude today. 

The conflict between theoretical physics and the ob- 
servational data can be alleviated by working within the 
framework of dynamical dark energy [2T| - [23| . This afore- 
said idea has led to a wide variety of dark energy mod- 
els such as quintessence d^-HII , k-essence f32|- [35[l. g uin- 
tom [36l - l45| , and holography dark energy (HDElpGl-liF 
In particular, the latter model was discussed extensively 
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during the last five years IM|59|, |59|, |6lH63l . 

The HDE model has its physical origin in the holo- 
graphic principle as well as some features related with 
string and quantum gravity theories [6J-|68| . The under- 
lying postulate can be stated as follows [68| : the number 
of degrees of freedom in a bounded system should be finite 
and is related to the area of its boundary. This principle 
also suggests that the ultraviolet (UV) cutoff scale of a 
system is connected to its infrared (IR) cutoff scale. In 
the case of a system with size L (IR length) and ultra- 
violet cutoff A without decaying into a black hole, it is 
required that the total energy in the region of size L 
should not exceed the mass of the black hole with the 
same size, thus, L^Pj^ < LM^p being Mp the reduced 
Planck mass whereas the UV cutoff scale is defined as 
A = py^lill. The largest L allowed is the one which sat- 
urates the above inequality and leads to an holographic 
dark energy given by pA = ic^ M^pL~^ , where c is a nu- 
merical factor. Hence, this principle connects the dark 
energy based on the quantum zero-point energy density 
caused by a short distance cutoff A with an IR cutoff ^ggl 
that is usually taken as the lar ge s cale of the universe, for 
instance, Hubble horizon j46l l47j. particle horizon [47| . 
event horizon [U or generalized IR cutoff [69l - l80| . 

A natural arena for investigating the coincidence prob- 
lem is consider a phenomcnological a ppr oach where dark 
matter interacts with dark energy [SlHSlI- From the 
observational point of view, an interacting dark sector 
is completely compatible with the current observations 
of standard candles and WAMP data [11,191. In the 
present paper, we show how it is possible to get a phys- 
ically viable model based on a new holographic dark 
energy density that interacts with dark matter. More 
precisely, it turns to be dark matter pc feels the pres- 
ence of dark energy px through the gravitational ex- 
pansion of the universe plus an exchange of energy be- 
tween themselves. Based on the holographic principle, 
we propose a dark energy model where the quantum 
zero point energy density pA is equal to the dark en- 
ergy density px being L an IR cutoff that will be re- 
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lated with a cosmological length. As a result of this, 
we take = p/^ = 3c^A/^pL^^ = 3c'^APpf{p,p) where 
fip,p) is an arbitrary positive function. This gives rise to 
self-interacting holographic dark energy models (SIHDE) 
where oc f{p,p), indicating that there is a coupling to 
the dark matter component. The new holographic dark 
energy model assumes a generalized IR cutoff L that de- 
pends on the total dark sector density p — Pc + Px and 
the pressure of the mixture p ~ Pc + Px- 

Several works have been devoted to obtain cosmologi- 
cal constraints in the case of Ricci scalar cutoff [tsI - ItH or 
generalized versions of this one [tIjIt^IsOI- For example, 
the joint analysis of the 307 union sample of SNIa, to- 
gether with CMB shift parameter given by WMAP5, and 
the BAO measurement from SDSS, suggest that the holo- 
graphic Ricci dark energy exhibits a quintom-like phase, 
so it leads to a new model consistent with the current ob- 
servation because the equation of state for the Ricci dark 
energy can cross the phantom line [73j . Using the gen- 
eral framework presented in jo^ , which suitably describes 
and unifies the dark sector with an exchange of energy, 
we will investigate a cosmological scenario where dark 
interacts with SIHDE. After that, we will confront our 
results with the current observational data and compare 
with the ACDM model. In the last section, we summarize 
our main results and conclude. 



II. EVOLUTION OF THE DARK 
COMPONENTS 

We consider a flat FRW universe filled with two com- 
ponents, dark matter and SIHDE with energy densities 
Pc and Px, respectively. We also assume that the equa- 
tions of states arc lOc = Pc/ Pc and = Px/Px, whereas 
the Einstein equations read 



Px, 



(1) 



-apc - /3px, 



(5) 



and comparing Eqs. ^ with Eq. ([5]), we obtain a com- 
patibility relation 



{a-ujc-l)r + 13 ■ 



(6) 



between the equation of state of both components and 
its ratio r = pd Px- This relation allows us to use the 
Eq. ([5]) with constant coefficients a and j3 instead of 
the Eq. ^ with non-constant coefficients. After solving 
the linear system of equations @ and ([5]) , we obtain the 
energy density of each dark component as functions of p 
and p' 



fip + p' 



ap + p 



(7) 



where A = a — /3 is the determinant of the linear system 
of equations. At this point, we introduce the interaction 
term, 3HQi, between the dark components by splitting 
the Eq. ([5]) in the following way 



apc 



-Qi 



I 

Px 



PPx = Qi. 



(8) 



After differentiating the first Eq. ([7]) and combining with 
Eq. we find a second order differential equation for 
the total energy density: 



p" + {a + /3)p' + a^p = QiA. 



(9) 



Once the interaction term Qi is selected and replaced 
in ([9|), the total energy density p of the dark sector is 
determined by solving the source equation 0. Having 
obtained p, we are in position to get pc and px from Eq. 
([7]). calculate the scale factor by integrating the Fried- 
mann equation ([1]), and find the equation of state of the 
mixture from the relation p = —p' — p. In the case of 
pressureless dark matter (cjc = 0), the equation of state 
of dark energy ^ is given by 



P'c + P'x + i^c + l)/5c + {(^x + i)Px = 0. 



(2) 



Here H = a/a stands for the Hubble expansion rate, a 
is the scale factor, and ' means derivative with respect 
to the variable rj = ln(a/ao)^ being ao the scale factor 
today. From Eqs. ([1]) and ^ the total pressure becomes, 
p = —p' — p, hence the SIHDE, p^ = f{p,p) turns p^ = 
f{p,p'). As already mentioned in the introduction, p^ is 
related with the UV cutoff, while L ~ f^^^^ is related 
to the IR cutoff. We now consider the simplest case of a 
linear SIHDE, 
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Px = 



a- /3 



(P + ap) , 



(3) 



where a and /3 arc both free constants. Rewriting Eqs. 
(H]) and © as 



Px, 



(4) 



= (a-l)r+(/3-l), 
so it becomes linear in r. 



(10) 



III. INTERACTING HOLOGRAPHIC MODEL 

In the present section, we are going to examine a pro- 
posal where the interaction term Qi is a general linear 
combination of pc, Px, P, and p' |94| 

Qi = ci ^ p + C2{uJs + 1 - ajPc- 



( , 1 fl> (w, + 1 - a) + 1 - /3) 

C3{uJs + 1 - li)px - C4 -. — v-r p ■ 11) 

[ujs + 1)A 

Here lo^ is a free constant parameter and the coefficients 
Ci fulfill the following condition ci + C2 + C3 + Ci = 1 in 
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order to assure the existence of stable power law solutions 
a = t2/3("=+i)||]. The case with C2 = Cg = C4 = was 
examined in |57l| . |85| , l93. The case ci = C2 = C4 = 
was analyzed in jOSMOq , |99l4l0l| . The linear interaction 
Qi oc p', ci = C2 = C3 = 0, was introduced in [9J] and 
now it is considered here for its study. 

Using Eqs. (O wc can rewrite the interaction ((TT|) as 
a linear combination of p and p' only, 



Qi = 



up + {us + ly^ju -{oJs-a + l){oJs -f3 + l)]p' 

A 

(12) 

where the parameter u is defined in terms of Wg, a, and 
/3 as follows: 

u ^ ciiiVs - a + l)(ws - /3 + 1) 



-C2/3K-a + l)-C3a(cj, -/3 + 1). (13) 



Replacing the interaction term (|12|) into the source 
equation we obtain a linear differential equation 

p" + {u;, + l)-^[{ujs + lf + a(3-u]p' + iap-u)p = 0, (14) 

whose characteristic polynomial roots are 

_i_ f3a — u 



7 = Us + 1, 



UJs + 1 



(15) 



We restrict our analysis to the case with positives roots 
in order to avoid phantom dark energy, then we choose 
< + 1 < 7+. Solving the Eq. dH]), we obtain 
the total energy density in terms of the scale factor and 
consenqucntly the effective pressure: 

p = feifl-s-r^ + 62a-3("=+i\ (16) 

p = (7+ - l)6ia~37+ ^ (17) 

From d?]) and ([TB]) . we get the dark matter and dark 
energy densities as a function of the scale factor 

(7+ - /?)bia-3^^ + [LPs -13 + l)62a''^^"°+^) 

Pc - 7 , (loj 



Px 



{a - 7+)6ia-37+ + (a - _ 1)620"^^"=+^^ 



(19) 

At very early times, dark matter and dark energy den- 
sities ([T5)) - behave as a~^^^ with a constant ratio 



(7+ - - 7+) while p ~ (7+ - 1)& 



la 



-37+ 



However, at late times the effective fluid, dark matter, 
and dark energy have the same behavior with the scale 
factor, namely, p ~ pc ~ Px — bia~^'-'^''^'^\ leading to 
n ~ (ws - /3 + l)/(a -1-ujs) and p ~ a-3("=+i). The 
aforesaid facts indicate that the interaction term Qi is 
a good candidate to represent adequately an interacting 



dark sector because the ratio dark matter-dark energy 
r has enough parameters to adjust the cosmological ob- 
servations and it also alleviates the so called coincidence 
problem. It is important to emphasize that the mutual 
exchange of energy between the dark components makes 
that their usual behavior with scale factor change radi- 
cally; we distingush in the dark densities two terms a~^'^ 
and a~'^("=+^'. In fact, wc will consider the case with 
7^ = 1 in order to get pressureless dark matter at early 
times. 



IV. OBSERVATIONAL DATA ANALYSIS 

In this section, we will perform some qualitative cos- 
mological constraints for the SIHDE model interacting 
with dark matter through the interaction term Qi pro- 
posed in last section. In order to do that, we start by 
constr aining the parameter space with the Hubb le data 
H{z) [T07t.|109|. and SNe la observations ^0^. The 
H(z) test was proba bly fi rst used to constrain cosmo- 
logical parameters in lll| and then in a large number of 
articles [79l.[80l [TT2l - ll23t . The statistical metho d requires 
the compilation of the observed value -ffobs [l07j |. |l09| 
and the best value for the present time z = taken from 
[lOSj l . Table U shows i?obs at different redshift with its 
corresponding la uncertainty and the reference where 
this value was reported. 



z 


H{z) 


la 


reference 




km s"^ Mpc"^ 


uncertainty 




0.000 


73.8 


±2.4 


[108] 


0.090 


69 


±12 


[109] 


0.170 


83 


±8 


[109] 


0.179 


75 


±4 


[12^ 


0.199 


75 


±5 


\m 


0.270 


77 


±14 


[109] 


0.352 


83 


±14 


[124] 


0.400 


95 


±17 


[109] 


0.480 


97 


±62 


[107] 


0.593 


104 


±13 


[124] 


0.680 


92 


±8 


\m 


0.781 


105 


±12 


[12^ 


0.875 


125 


±17 


[124J 


0.880 


90 


±40 


[107] 


1.037 


154 


±20 


[124] 


1.300 


168 


±17 


[109] 


1.430 


177 


±18 


[109] 


1.530 


140 


±14 


[109] 


1.750 


202 


±40 


[1091 



TABLE I: Hubble data H^^^{ 



Zi 1 vs. rC' 



dshift . 



From Eqs. ([T|), ([H]), and ([H]), we can write the Hubble 
function in terms of the effective equation of state a; = 



4 



H{z) = i/o 



pn.^, = -2H/3H^ as follows 



(w, -a;o)(l + 2;)''+a;o(l + 2)'^ 



(20) 

where wq = "^^co + /^^^xo - 1, ^xo = Pxo/^Hq, flco = 
Pco ^^"^ their present values whereas the flatness con- 
dition today reads flco + ^^^o = 1- Taking into account 
the transition point Zacc, i-e. the moment where the uni- 
verse begins to accelerate or where the deceleration pa- 
rameter q — —a/aH^ vanishes, the Eq. ((20|) depends on 
{Hq,z^cc,^s) parameters: 



H{z) = 



\3/2 



v/-3a;, -f + (H-Zacc)"3'' 



(1 + ^) 



3(a;,-l-l) 



(1 + Zacc)3("^ + 1) 



(l + 3w,) 



(1 + Zacc)^ 



1/2 



(21) 



From Ea.([2T|) we see that the model has only three in- 
dependent parameters (ifo, ^acc, Ws) in order to be com- 
pletely specified. The remaining parameters a, /3, fico 
or n^o are included in the transition point Zacc through 
the equation of state wq. We now proceed in the fol- 
lowing way: we perform a statistical analysis on the 
(Hq, Zacc, ^s) parameters, confronting their best fit values 
with the recent available data and then, we will perform 
the same Hubble test using the expression ((20)) , to obtain 
constraints on the others parameters, a, /3, and fJcO- The 
second approach will give us the most favored SIHDE for 
the Hubblc's data. 

The probability distribution for the 0-parameters is 
P{e) = 7Ve-'^'('')/2 (see e.g. being M a normal- 

ization constant. The parameters of the model are esti- 
mated by minimizing the function of the Hubble data 
which is constructed as 



(t2(z,) 



(22) 



where 6 stands for the cosmological parameters, Hohs{zk) 
is the observational H{z) data at the redshift z^, a{zk) 
is the corresponding la uncertainty, and the summation 
is over the 19 observational H{z) data. The Hubble 
function is not integrated over and it is directly related 
with the properties of the dark energy, since its value 
comes from the cosmological observations. Using the ab- 
solute ages of passively evolving galaxies observed at dif- 
ferent redshifts, one obtains the differential ages dz/dt 
and the function H(z) can be measured through the re- 
lation H{z) = —{l + z)~'^dz/dt. The function reaches 
its minimum value at the best fit value 9c and the fit is 
good when Xmuiifc) /[^ — 1j where n is the number 
of parameters [l03| and N counts the observational data 
points that in our case correspond to 19 points. 

In the first approach, the parameters of the model arc 
9 = iHo,Zacc,oJs) therefore the 68.3% (Icr) or 95.4% 




FIG. 1: Three-dimensional C.L. assocaited with Icr, 2a for Hq^ u^, and 
2acc parameters. The point indicates the best fit observational values, 
namely, Hq = 71.839 kms~^ Mpc"^, z^cc = 0.7831 and ojj, = -0.950. 



0.415 



1.04 



-0.9.'' 




0.78 



FIG. 2: Two-dimensional C.L. assocaited with la, 2a for ixJs and Za.cc 
parameters, after the marginalization over over the parameter Hq was 
done. The point indicates the best fit observational value obtained with 
the H{z) function. 



{2a) confidence levels (C.L.) made with the random 
data fulfill the inequalities x^W - Xmin(^c) < 3.53 
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FIG. 3: Constraints in the cj^-ficO plane. Elliptical two-dimensional 
C.L. associated with \a and la error bars. The dot indicates the best 
fit observational value obtained with the H{z) function. The shaded 
zone is excluded because we try to avoid phantom dark energy in our 
model. 



103|-|l05|. Taking into account for the sta- 



tistical analysis, we obtain the best-fit values a = 1.15, 
P = 0.023, and fl^o = 0.243 with xLn = 14.2623 along 
with Xd o f = 0.95 < 1. Fig 12] shows two-dimensional 
C.L. in the ajg-i^co plane whereas the other parameters 
are taken as priors, namely, we fix a = 1.15, /? — 0.023, 
and Ho = 71.84 kms'^Mpc-^ Then, the best-fit val- 
ues together with their error bars arc iMs = —0.95^^^03 
and rjco = 0.244+"o°|^g. We would like to use the best- 




or x^i(^) - Xmin(^c) < 8.02, respectively. Fig H] shows 
the C.L. associated with la and 2a error bars in the 
Hq — Zacc ~ space; we find the best-fit values at Hq = 
71.839 kms-^Mpc-\ Zacc = 0.7831 and a;^ = -0.95 
corresponding to a Xmin = 14.26 along with Xdof — 
Xmin/(-^ — Ti) = 0.790 per degree of freedom. We remark 
that our estimations of the actual Hub ble p arameter 
agree with the median statistics made in [IQJ], namely, 
our value meets within the la interval obtained with the 
median statistics. Ho = 6 8 ± 5 kms~^Mpc~^, or with 
the analysis performed in |105| about the impact of Ho 
prior on the evidence for dark radiation. On the other 
hand, we obtain C.L. in the Zacc — plane obtained after 
having marginalized the joint probability P{Ho, .Jacci '^s) 
over Ho (see Fig 12]). As usual, in the case of two pa- 
rameters, 68.3%, 95.4% C.L. are made of random data 
sets that satisfy the inequality x^iO) — \l ,;„(Or) < 2.3, 
X^(^) - xLini^c) < 6.17, respectively [l^. The shaded 
band corresponding to Ws < — 1 is excluded in our model 
in order to avoid phantom dark energy. The constraint 
on the critical redshift is Zacc = ^■'^^^-037: such value are 
in agreement with zt = 0.691q^3 reported in [l25l |- 126 |. 
and meets within the 2 a C .L obtained with the super- 
novae (Union 2) data in [l26j | . The critical redshift Zacc = 



0.78+Vf7 is also consistent with Union 2 + BAG + CMB 
data jl06l |. For the other parameter the statistical anal- 
ysis leads to ujs = -0.95"t°i^27- 

In order to get some physically relevant bounds on a, 
/3, and fico parameters, we now use the expression (|20p 
and take as prior Ho = 71.84 kms~^ Mpc"-'^ which is in 
agreement with the median statistical constraints found 



FIG. 4: The plot of the theoretical distance modulus (solid line) v ersus 
the redshift. The observational data (point) was taken from Riess fTol 
and covers two different regions of redshift, z < 0.1 and 0.1 < z < 1.5. 
We obtain that the best-fit values, obtained for the interacting model, 
arc in agreement with the supcrnovac data. 



fit values Ws = —0.95 and fico = 0.244 to calculate the 
magnitude redshift relation for standard candles and con- 
trast with the supernova data. As it is well known the 
observations of SNe la have predicted and confirmed that 
our universe is currently passing through an accelerated 
phase of expansion. Since then, the observational data 
coming from these standard candles have been taken very 
seriously. It is commonly believed that by measuring 
both their redshifts and apparent peak flux gives a direct 
measurement of their luminosity distances and thus SNe 
la provide the strongest constraint on the cosmological 
parameters. The theoretical distance modulus is dcflned 
as 



/x(z) 5 logioPi + ^0, 



(23) 



where no = 43.028, and 2?l is the Hubble-free luminosity 
distance, which for a spatially flat universe can be recast 
as 



Vl{z) = (l + z)Ho f 
Jo 



dz' 
H{?)- 



(24) 



Replacing the best-fit values of Ho, lOs, and fico in 
Eas. (pni) - (|M)) wc get the theoretical distance modulus 
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/j,{z) for our model (see Fig. |4]) whereas the observa- 
tional data with their error bars, ^obs{zi), are taken from 
[l08l |. As we can see from Fig HI our model exhibit an 
excellent agreement with the observational data, at least 
in the zones corresponding to small redshifts [z < 0.1] 
and large redshifts [0.1 < z < 1.5]. 



A. The age problem 

We now turn our attention to the age problem, namely, 
the u niverse cannot be younger than its constituents (see 
127|). For example, the matter-dominated FRW uni- 



verse can be ruled out because its age is smaller than 
the ages inferred from old globular clusters. The age 
problem becomes even more serious when we consider 
the age of the universe at high redshift. Now, there are 
some old high redshift objects (OHROs) discovered, for 
instance, the 3 .5 Gyr ol d galaxy LBDS 53W091 at red- 
shift z = 1.55 [miili, the 4.0 Gyr old galaxy LBDS 
53W069 at redshift z = 1.43 [Tsol. the 4.0 Gyr old ra- 
dio galaxy 3C 65 at z = 1.175 |13l| . and the high red- 
shift quasar B1422-I-231 at z = 3.62 whose best -fit age 
is 1.5 Gyr with a lower bound of 1.3 Gyr jl32| . Also 
the old quasar APM 08279-^5255 at z = 3.91, whose age 
is estimated to be 2.0 — 3.0 Gyr [l33l Il34j . is used ex- 
tensively. To assure the robustness of our analysis, we 
use the most conservative lower age estimate 2 .0 Gyr fo r 
the old quasar APM 08279-f5255 at z = 3.91 [1331 . [T33 |. 
and the lower age estimate 1.3 Gyr for the high redshift 
quasar B1422-f 231 at z = 3.62 |l32{ . Many authors have 
examined the age problem within the framework of the 
dark energy models, see e.g. |l27l |. [l35]- 141, 142], and 
references therein. The age problem within the context 
of h olographi c dark energy model was explored in jl4C 



and In this section, we would like to consider 



the age problem for the SIHDE model with linear inter- 
action. 

The age of our universe at redshift z can b e ob tained 
from the dimcnsionless age parameter ( (l27j ] , jl36j ] ) 



r°° dz' 
Uz) = Hotiz) = Hoj^ (7TTWM 



(25) 



At any redshift, the age of our universe should be larger 
or equal than the age of the old high redshift objects 



Tziz) > Tobj = Hatobj, or S{2 



Tobj 



> 1 



(26) 



where tobj is the age of the OHRO. It is worth noting 
that from Eq. (|25|) . Tz{z) is independent of the Hubble 
constant Hq. On the other hand, from Eq. ([25)) . Tobj is 
proportional to the Hubble constant Hq that we consider 
as Ha = 71.84 kms-^Mpc-^. 

In Table HI we show the ratio S{z) = Tz{z)/Tobj 
at z = 3.91,3.62,1.55,1.43,1.175 taking to account the 
best-fit values obtained in the last section. We obtain 
that r,(z) > Tobj{z) at z = 3.62,1.55,1.43,1.175 but 



Tz{z) < Tob]{z) at z = 3.91, so the old quasar APM 
08279-1-5255 cannot be accommodated as the others old 
objects. Perhaps, the age crisis at high redshift in the 
case of dark energy holographic models jl40j ]. [? ] could 
be alleviated by taking into account another type of in- 
teraction. This fact will be explored in a future research. 



5(3.91) 5(3.62) 5(1.55) 5(1.43) 5(1.175) 



0.854555 1.19781 1.20467 1.12826 1.31645 



TABLE II: It shows the ratio S(z) = T^(z)/Tabj at z = 
3.91,3.62,1.55,1.43,1.175 for the best-fit values obtained with the 
Hubble data. 



B. Kinematic analysis 

FigO shows the behavior of the deceleration param- 
eter with redshifts. Using the values ilco = 0.244 and 
^xo = 0.75, we obtain that the deceleration parame- 
ter vanishes at Zacc = 0.78, so the universe enters the 
accelerated phase earlier than the ACDM model. Re- 
garding the effective equation of state, it stays in the 
range —1 < uj{z) < for z > 0, more precisely, oj{z) 
starts as non relativistic cold matter, decreases rapidly 
around z = 2 and then ends with the asymptotic value 
LJs = —0.95. The dark energy equation of state stays 
in the range —1 < Wj-(z) < also. The ratio dark matter- 
dark enery, 

^ _ ujoiP - CO, - 1) + (1 - /3)(^o - c^,)(l + zf^' 

wo(ws + l-a) + (a-l)(a;o-u;s)(l-hz)3'^=' ^ ' 

evaluated at the best fit values a = 1.15, /3 = 0.024, 
LOs = —0.95 and ujq = 0.78 indicates that interaction Qi 
helped to alleviate the coincidence problem. 



V. SUMMARY AND CONCLUSIONS 

In this paper we have considered a flat FRW universe 
composed of an interacting dark matter and SIHDE. We 
have shown that the compatibility between SIHDE and 
the conservation equation gives a constraint between the 
equations of state of the dark components. We have se- 
lected linear SIHDE and a linear interaction in the dark 
sector and find that this model describes properly the 
evolution of both dark components. We have also shown 
that a general linear interaction, Q;, is a good candidate 
for alleviating the cosmic coincidence problem. 

Taking into account the Hubble data (see Table |l| and 
using the statistical method, we have obtained the 
best-fit values at Hq = 71.839 kms~^Mpc~^, Zacc = 
0.7831 and = -0.95 along with xlo.f = 0-790 < 1 
per degree of freedom (see Fig. [T]). The value of 
Hq is i n ag reement with the one reported in the lit- 
erature 11081 or with the median statistical constraints 
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0.0 0.5 1.0 1.5 2.0 
z 



FIG. 5: Plot of deceleration parameter q{z) taking into account the 
best-fit values fl^ = 0.75, a = 1.15, P = 0.024, and lu^ = -0.95 
(magenta, solid line). It also shows the deceleration parameter for the 
AC DAI model (green, dashed line). 



(Union 2) data in jl26j |. I t is also consistent with 
Union 2 + B AO + CMB data . Using as priors a = 
1.15, /3 = 0.023, and Hq = 71.84 kms'iMpc-\ we 
build two-dimensional C.L. in the iOg-^lcO plane and es- 
timated the best-fit values ujs = — 0.95^;i^^Q3 and flco = 
0.244+_°o°l89 (seeFigEl). 

Replacing the best- fit values of Hq, uis, and rJcO; we 
obtained the theoretical distance modulus iJ,{z) for our 
model (see Fig. |4]) a nd co nfronted with supernovae data 
IJ-obsizi) taken from jlOSj . Fig. S] shows that our model 
exhibit an excellent agreement with the observational 
data. Besides, we have found that the age crisis at high 
redshift cannot be alleviated because the old quasar APM 
082794-5255 at z = 3.91 [l33, 134j seems to be older than 
the universe; so it will be needed to consider other kind 
of interaction (cf. Table |ll| or perhaps to propose a non 
linear SIHDE for exploring this issue. Finally, we have 
found that our model enters the accelerated phase faster 
than the ACDM model (see FiglS]). Concerning the ef- 
fective equation of state and the dark energy equation 
of state, we found that both do not cross the phantom 
divide line (4l| . In a future research, we are going to ex- 
plore the linear SIHDE proposal where the dark sector 
is also coupled to a radiation or baryonic term; we will 
examine the changes i ntro duced in the behavior of dark 
energy at early times jl46l | . 



found in jl04| . [l05l |. Having marginalized the joint 
probability P{Ho, z^cc^s) over Hq (see Fig 12]) we build 
two dimensional C.L. and obtained the best-fit values 
with their Icr error bars, namely, Zacc = 0.78~'^q^3-7 and 
ujs = — 0.95^i^27- The critical redshift is in agree- 
ment with zt = 0.69to i3 reported in [l2i]-[l2i], and 
meets within the 2cr C.L obtained with the supernovae 
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